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Abstract. 

We investigate the linear stability of topological black hole solutions to four-dimensional su(AT) Einstein-Yang- 
Mills theory with a negative cosmological constant. We here extend recent results in the field which prove the 
existence of hairy black hole solutions to such equations [l], and the stability of their spherically symmetric 
analogues [2]. We find the analysis in [2] carries over very similarly, with some important differences in the final 
stages. Nevertheless, we establish the existence of non-trivial solutions stable under linear perturbations, in a 
sufficiently small neighbourhood of some existing trivial solutions; in fact, stable topological solutions turn out to 
be likely more abundant in the parameter space than their spherically symmetric analogues. 
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1. Introduction 

Recent work has shed much light on hairy black holes and solitons - that is, solutions to four-dimensional Einstein- 
Yang-Mills (EYM) theory. From the first solutions discovered in four-dimensional su(2) asymptotically flat space 
by Bartnik and McKinnon [3j in the soliton case and Bizon [3] in the black hole case, interest has been high 
in investigating many generalisations of their work, for instance: to spaces of higher dimension or alternate 
topologies EH2; to ‘dyonic’ solutions, characterised by a non-zero electric gauge sector [HlfTPl (which only exist 
in spaces where A < 0); and to extensions of the gauge group to su(N) or to abitrary compact gauge group [7X1172] . 

Of central importance to the physicality of these solutions is the question of their stability. We may prove 
that certain solutions exist, but if we are to prove their longevity in the universe, they must be stable under linear 
perturbations at the very least; otherwise we cannot be sure that if disturbed, the black hole hair will not radiate 
away or collapse inwards. 

The first cases which were considered were all for asymptotically flat space. However, it was proven by 
Straumann and Zhou in [73 that the original solutions discovered in su(2) were unstable under linear spherically 
symmetric perturbations. It was proven that the number of unstable modes of the n th solution (indexed by n, the 
number of zeroes (‘nodes’) that the gauge field function to possesses) is 2n, with n in each of the two decoupled 
‘sectors’ (see later for details) and that for A = 0, the gauge field must possess at least one zero [71H7B] . Zhou 
performed non-linear analysis on solitons in su(2) and again showed them to be unstable, in that either hair radiates 
to infinity or else collapses to produce a Schwarzschild black hole m- Finally, Brodbeck and Straumann proved 
instability for arbitrary gauge group in the asymptotically flat case [18] . For a comprehensive review of stability in 
these cases, see [751 . 

It also makes physical sense that stability is directly related to the sign of the cosmological constant A 
specifically, that solutions are only expected to be stable for A < 0. This is a consequence of the ‘balancing act’ 
between the repulsive forces set up by the gauge field and the attractive force of gravity; in the open geometry 
associated with A > 0, hair is free to radiate away to infinity if disturbed. However, A < 0 corresponds to a closed 
geometry, meaning that the black hole hair is ‘trapped’. This is reflected in the nature of the solution space in each 
case: it is known that, in general, asymptotically flat hairy solutions occur only at discrete values of the parameters 
concerned, in other words their solution set does not form a continuum [20H23| . whereas for asymptotically anti-de 
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Sitter (adS) space, solutions occur in continuous ranges of the parameter spaces and so are more plentiful | [2iU2T] . 
This gives us hope that small perturbations in existing solutions will simply yield new solutions ‘not too far’ from 
the original, and hence these solutions may be stable. 

Many solutions have been found for A < 0 in a variety of cases: see [2S] for a recent review. The stability of 
some su(2) solutions in asymptotically adS space is proven for |A| —> oo in [2M50] . In addition, recent work [2] shows 
stability for some solutions in spherically symmetric su(.ZV) EYM theory. Given that we have recently established 
the existence of their topological analogues [I], it seemed natural to extend the work to cover a stability analysis 
and thus fill a gap in the research. Upon examination, it transpires that the analysis of the topological case differs 
from the spherical case in some small but vital details, which bear future consideration. 

The purpose of this work is to examine those details, and so because of this we direct the reader towards [2] 
and note that since the basic analysis is so similar, we shall only concentrate here on the cases for which k ^ 1 
and focus solely on where the cases differ. We emphasise, this work is written to be read alongside [2] as we rely 
heavily upon it. The outline of this paper is as follows. We begin by outlining all relevant formulae and ansatze. 
We briefly review the static field equations, boundary conditions and trivial solutions, and move quickly on to the 
problem of stability. We find again that our choice to consider purely magnetic solutions allows the solution space 
to decouple into two sectors, which we take in turn - we call these the ‘sphaleronic’ and ‘gravitational’ sectors. 
We find that the topology of the manifold makes little difference until the end of the analysis, when it becomes 
important: nevertheless we are able to use a theorem of Arnann and Quittner ’ S m to prove the linear stability 
of non-trivial solutions in a small neighbourhood of existing trivial solutions, and in addition to this, we find that 
stable solutions appear to be more plentiful in the strictly topological case compared to the spherically symmetric 
case. 

2. Topological su(N) Einstein-Yang-Mills theory 

In this section we give the general mathematical background we shall need to describe and model topological black 
holes. By way of review, we will also show the equilibrium field equations, boundary conditions and embedded 
solutions. 

2.1. Ansatze 

The action used for the four-dimensional su(N) Einstein-Yang-Mills (EYM) theory, with the cosmological constant 
A < 0, is 

Seym = \ j d A x^[R - 2A - Tr F^F^], (1) 

where R is the Ricci scalar, is the non-Abelian gauge field and Tr signifies the Lie algebra trace. Throughout 
the paper the metric has signature (— , +, +, +) and we use units in which 47rG = 1 = c. 

Varying the action gives the field equations 

2 = R^v — -g^R + A g^, 

0 = V x F^ + [A x ,F^}, 

where the Yang-Mills (YM) stress-energy tensor T and the field strength tensor F^ are given by 

Tfiu = Tr [F^ X F X - \g^F Xa F Xa ] , F^ = d^A v - d v A M + i[A^, A„], (3) 

and we have rescaled so that the gauge coupling constant g = 1. Note also that we have employed the usual Einstein 
summation convention where it is understood that summation occurs over repeated indices. 

In this paper we are interested in static, topological black hole solutions of the field equations ©, specifically 
for spaces regularly foliated by 2D (spacelike) hypersurfaces of constant and unit- or zero-magnitude Gaussian 
curvature fc, and hence (following [T]) we write the metric in standard Schwarzschild co-ordinates as 

ds 2 = -gS 2 dt 2 + p~ 1 dr 2 + r 2 d9 2 + r 2 fl(Q)dtf , (4) 

where p and S depend on t and r alone. For convenience, we may take 

2m (t, 7') A r 2 

3~ 


g(t, r ) = k 


r 


(5) 
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The function fk{0) is used to endow the metric with the correct topology in the angular portion, and depends on 
our value of k such that 

{ sin# for k = 1, 

9 for k = 0, (6) 

sinh # for k = — 1. 

It should stressed that we will solely be concerned with the non-spherically symmetric cases, i.e. k = 0,-1. In pQ, 
we used Wang’s theorem [52] to construct the following topological gauge potential: 


A = Adt + Bdr + -(C - C + )d# - - 


(<C + C*)f k (0)+D ^ 


d(f>. 


(7) 


In the above A, B , D are all (TV x TV) traceless, diagonal matrices, with A and B imaginary, D real; C is a complex 
(TV x TV) upper-triangular matrix with non-zero entries only immediately above the main diagonal; and C f is the 
Hermitian conjugate of C. (See jl] for details.) In comparison with (2.8]3, the Einstein tensor G^ that we derive 
from the metric ansatz has components given bj@: 

G tt = (fJ-'r - k + n ), G tr = -—, G rr = l 2 (p'Sr + 2 S’pr - kS + pS). (8) 

r z pr [iSr^ 

It is unnecessary to state the angular components of the Einstein equations since if © hold, then the other 
components do automatically due to the Bianchi identities. We also need not state the components of since 
three of the components of the field strength tensor ( F tr , F t g and F r g) are the same as in (2.10) of [2j, and three 
are slightly different: in F t( j >1 F r< ^ and Fg the ‘sin#’ is replaced by /&(#). 

Finally, we comment on the equilibrium field equations, the boundary conditions, and the embedded solutions 
here. The static field equations and symmetries are the same as (2.17 - 2.22) in )2] but for the following changes: 

1 N 

• The function fl (2.18) is altered to n = —j ^ (ipj 

r 3 =i 

• The function Wj ( 2 . 20 ) becomes W 3 p/. 

Wj, k = k — Wj + - (uj +1 + Wj-i) ■ ( 9 ) 


-^ 2 _ 1 -fc(TV + l-2j)) 2 ); 


The stated boundary conditions (2.23 - 2.26) are identical but for the alteration of (2.24) to m{i'h) = 


M 

6 


There will be some changes in the higher order parameters which involve k but in this paper, as in [2j, we are only 
interested in the leading order behaviour of the field variables. As for the embedded solutions, we note that the 
SAdS solution (2.27 - 2.28) is unavailable for k ^ 1, but the other two solutions remain: the Reissner-Nordstrom- 
Topological-anti-de-Sitter (RNTadS) solution (2.29 - 2.30), which has an altered magnetic charge (2.31) given 
by 


fc 2 TV(TV 2 - 1) 
Qm =- 2 -> 


( 10 ) 


and the embedded su(2) solution (2.32), defined by letting u>j = \/j(N — j)w. The only other detail which it may 
be important to state is that the bound given by requiring m'(vh) > 0 leads to a minimum event horizon radius for 
k = —1, given by 

Vh\ > a/|A|(2 m'{r h ) + 1), i.e. \r h \ > (|A| - 2r 2 n(r /l )) 2 , (11) 

which in turn implies a minimum bound on |A|, generalising a result in [5]: 


|A|>l(l + 2r 2 n (r h )). (12) 

r h 

J For brevity, all section numbers etc. which are in bold will be understood to refer to appropriate sections of [2]; equations from [2] 
will be simply given in the form (X.XX). 

§ Note that throughout this work, we denote ‘ to mean d/dt and ' to mean d/dr. 
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We wish to examine the behaviour of solutions to the field equations under linear perturbations, so all Held variables 
(/x, to, S , cjj, ctj , /3j, 7 j) are rewritten the form 

= f{r) +Sf{t,r) (13) 

(where from now on, a bar “ represents the equilibrium solution variables, which are dependent on r alone, and 
where dy = qy- = 0), and the resulting perturbation equations are obtained. It should be noted that in our case, 
the perturbations will of course not be spherically symmetric, but will share the symmetry (planar or hyperbolic) 
of the angular portion of the topological metric (g]|. We are only interested in purely magnetic solutions; therefore 
we have a zero electric gauge field and we set the perturbations day = 0. In this regime, the equations decouple 
into two sectors: the sphaleronic sec<o?[jJ, containing only the YM perturbations 5f3j and S'jj; and the gravitational 
sector , containing only the metric perturbations Sp and SS and the YM perturbations Suij. Thus, these two sectors 
are considered separately. 

A final point to make before we do this is that Sections 2.3, 3.1 - 3.4, 3.6, 4.1 and 4.2 carry across almost 
identically, with two main (and small) differences, both of which are the same as the alterations in the field equations 
(see the list following equation (|8jl); the result of which is that only terms involving Wj t k actually make a difference 
to the final analysis. As we will show though this does make a fairly significant difference to the results. 


3.1. Sphaleronic sector 


The basic idea here was to use co-ordinate transforms to organise the variables S/3j and 5"/j into a single vector 
and thus to state their second order differential equations in a self-adjoint hyperbolic form (3.19): 

U^_ = -4/ (14) 

Then, for time-periodic perturbations, i.e. for which &_(t,r) = Re (e z<T, ^(r)), the perturbations will take the form 


= <t 2 \L (15) 

so that if U > 0 then the eigenvalues er 2 > 0, so that a € R, making e lat }&_(r) a time-periodic perturbation and thus 
stable. (If U = 0, we get no dynamics and the only perturbations correspond to gauge transformations - hence we 
get trivial stability.) 

We arranged the operator U into the form IA = x'X + V — G^h 2 G, where x and V are matrices in block-matrix 
form, G is the operator representing the Gauss constraint, and h = S y' . Due to the forms of the various terms 
involved, and the fact that the last term will vanish when applied to physical perturbations that satisfy the strong 
Gauss constraint, we ensured that U was hyperbolic, symmetric and real; and also positive, given only one constraint 
appearing in the matrix V: that is, 

W < 0, (16) 

for the matrix W (3.15). Following the analysis, this means the argument is structurally identical in our case; and 
in fact the only difference that survives the process is in W, and hence the matrix element V 22 (3.39). Therefore, 
the main difference in the analysis of Sections 2 - 4 is in the stability condition (fl6l) . dealt with in Sections 3.5 and 
4.3 ■ the so-called special cases. What we need to consider first is what (HI implies for those cases. 


3.1.1. Conditions for no sphaleronic sector instabilities Comparing our work with Section 3.1 - 3.4, the only 
difference is in the matrix W (3.15), where in our case 

W = / l 2 diag(W 1 , fe ,W 2 ,fc,...) J (17) 

for which Wjj. is the equilibrium part of Wy/c, defined in (0. Therefore for the matrix W to be non-positive 
we require Wj t k < 0, which gives us a set of inequalities that must be satisfied by the equilibrium solutions, for 
j = 1, ..., N —1 - these are the analogues of (3.42): 




k + 0 + • 


(18) 


|| This sector is so named because the original su(2) EYM solitons [3] and black holes [3] in asymptotically flat space possess instabilities 
in this sector [H11IS] analogous to the unstable mode of the Yang-Mills-Higgs sphaleron m- 
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The first thing to notice is that if we take k = 0,-1, the right-hand side of (1T51) is lower than if k = l, so naively 
we may think that stable solutions in these cases might be more abundant. Once again though, the challenge is to 
prove that we can find equilibrium solutions which conform to (1181) for all r. 

3.1.2. Special cases As in [2], we must next consider the embedded solutions discussed in Sections 3.5.2, 3.5.3. 

The RNTadS solution Here we let CCj = 0 for j = — 1. We find that the perturbations reduce down 

considerably, to the point where the only allowed solution for bound state perturbations is the trivial solution 
J/ = 0. Hence, there are no dynamics in the sphaleronic sector in this case, and the only bound state perturbations 
are pure gauge transformations. So the RNTadS solution is trivially stable here. This is identical to the k = 1 case. 


su(2) embedded solutions Here we let uoj = \/j(N — j ) for j = 1, ...,N — 1. In that case, 

W = h 2 (k-Cd 2 )T N - 1, (19) 

which means that the condition W < 0 implies 

Co 2 > k (20) 

for all r. Thus, all topological su(2) solutions immediately satisfy (1201) for all functions Co. 

Finally, we note that since these embedded su(2) solutions satisfy dl8l) . the analysis of Section 3.6 in [ 2 ] carries 
straight over (using instead the analogous existence propositions from [I to construct the nearby solution); and 
hence we can still find non-trivial solutions which are sufficiently close to existing embedded solutions, within some 
small enough neighbourhood of the parameter space, and which possess no instabilities in the sphaleronic sector. 

3.2. Gravitational sector 

The situation is somewhat similar to the sphaleronic sector, though the analysis is far more resistant to simplification 
in the general case. We begin by defining a vector do; = (daq, St 02 , ...dw,/v_i) T . After manipulating the perturbation 
equations to make the problem more tractable, we showed that the gravitational sector is naturally self-adjoint. 
Thus we expressed the YM perturbations SiOj in the following matrix form (4.8): 

— SCo = — d 2 t 6ui + AdSiv (21) 

Therefore, for all time-periodic perturbations such that Soj_(t,r) = e lcrt Suo_(r) , we have 

ct 2 <5 u > = -^Su+MSu. (22) 

Since (1221) is a a Schrodinger-like equation, the non-negativity of Ad implies that a £ R. and hence there are no 
unstable modes in the gravitational sector; so we wish to prove in each case that Ad > 0. 

As in the k = 1 case, the matrix A4 (4.9) can be separated into entries that are on the main diagonal (Adjj), 
entries just above and below the main diagonal (A4jj + 1 = Aij+ij), and all other entries (Adj t k), as follows: 


7j 02 A OC 

M jt j = - - 2u, 2 ) - —~Q(d r ,Oo 3 ) 2 - —(dr.CodW^Co, 


Ad 


jj +1 


fj,Sr 

— 2 ~U}jC 0 j + 1 - -r^-Q(dr„Uj)(dr,U)j+l) - J (Wj : kCjjdr,LOj +1 + Wj+l,kCbj + ldr,U>j) 

V /AID T T 


4 8 S - 

flo T T 


(23) 


where r» is the ‘tortoise’ co-ordinate introduced in [2] (and e.g. JM]), i.e. 


dr* 1 
dr fiS 


(24) 
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3.2.1. The RNTadS solution The first embedded solution we shall test is the RNTadS solution, which is found by 
letting Cbj = 0 for all j. Then we find that W 3t k = k , and the matrix Ad in (1231) reduces right down to (in analogy 
with (4.15)) 

n s 2 

M = -k^-1. (25) 

Let’s take the 2 cases of k which concern us: 

k = 0: Here, Ad is the zero matrix, and so the equation reduces to a 2 6uj_ = —d 2 6uj_. Since the right-hand side of 
this is still positive, then a 2 > 0, giving a e M; and therefore the RNTadS solution has a stable gravitational sector 
here. 

k = — 1: Now, the matrix Ad has positive eigenvalues, again meaning the right-hand side is positive giving cr € K, 
and therefore once more the RNTAdS solution has a stable gravitational sector. 

We note that both of these results are directly the opposite to the spherically symmetric case, and so we have 
result of some interest. It appears that we have a solution in which uij = 0 for all j, which as we have explained, 
we might not necessarily expect stability from - however, the RNTadS solution is stable for k ^ 1. We will return 
to this point in Section |4] 

3.2.2. su(2) embedded solutions Embedding su(N) in su(2) means making the following identification: 

=uj(r)y/j(N - j). (26) 

We begin by rewriting the non-zero elements of the original Ad matrix using this identification. We note in this 
case that Wjj.- = k — to 2 . Then we find we can write the matrix Ad (1231) more conveniently by separating the matrix 
into three separate terms (4.17 - 4.18): 

Ad = A/). + A /2 + A/ 3 . (27) 

However, A /2 and A /3 are identical to those in (2], thus neither contain any reference to k and so their positivity is 
assured; unconditionally for A/ 2 , and for A /3 under the condition that |A| —> 00. Therefore we only need to examine 

A/i: 

A/i = ^(w 2 -fc)I. (28) 

We bear in mind that for stability, we want this part to be non-negative, in which case we have 

w 2 > k. (29) 

This is the same as the condition for stability in the sphaleronic sector; and it is guaranteed to be fulfilled for 
k = 0,-1. Therefore A/"i > 0, Ad is non-negative, and thus all su(2) embedded solutions possess no unstable nodes 
in the gravitational sector provided |A| —> 00. Comparing the k = 1 case then, we see that this situation is similar 
but without any conditions on ujj , meaning our region of stable solutions in the parameter space should in principle 
be larger than in the spherically symmetric space. Note that the regime of |A| large was also used in [5] to prove 
stability in the gravitational sector for topological su( 2 ) solutions. 

We briefly summarise our results so far. In the sphaleronic sector, the RNTadS solution has no dynamics; 
the su( 2 ) embedded solution is unconditionally stable; and we can find nearby non-trivial solutions stable in the 
sector, i.e. they also fulfil the condition (1181) . Therefore, we only need to be concerned with the gravitational 
sector; in which the RNTadS solution is unconditionally stable and the su(2) embedded solution is stable for large 
|A|. All there is left to do is prove that non-trivial solutions exist which are stable in the gravitational sector, 
in a neighborhood of any su(2) embedded solutions. The rest of Section 4 employs a theorem by Amann and 
Quittner [31| that we will refer to as the Nodal Theorem. 

3.2.3. The Nodal Theorem This theorem (quoted in Section 4.4.1) concerns the eigenvalue problem for a certain 
Schrodinger-like linear differential operator, and determines the number of bound states of such an operator in 
terms of the nodes of the solution vectors of an auxiliary initial value problem (4.26). We note that Section 4.4 
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carries over virtually unchanged (ignoring Section 4.4.3 on solitons): in particular, equations (4.28) to (4.35) in 
Section (4.4.2) are identical to our case but for one occurrence of k in A ij.j (4.31), which becomes 


MjJ ~ 3 


k 3cUyoo “b 2 l,c 


^j+ l,oo) 


AO^ 1 ), 


(30) 


though this doesn’t change the leading order of the expression and hence the application of the Nodal Theorem is 
unaffected. Thus, it can be proven that su(2) embedded solutions possess no unstable modes in the gravitational 
sector. 

Finally it remains to prove that we can find non-trivial solutions in some neighbourhood of these embedded 
solutions that still satisfy the Nodal Theorem. This is dealt with in Section 4.5. The basic idea is that for the su(2) 
embedded solutions we prove that the Nodal Theorem is fulfilled as long as the function defined by 


5{p) = detll c (p) (31) 

has no zeroes over a particular range, where U c = \u\,U2, ■■■iM.n-i] (4.25) is the matrix of solution vectors to the 
eigenvalue problem (4.26), and the vectors Uj represent the gauge field perturbations Suij. We then must apply the 
theorem to prove that we can find non-trivial solutions nearby the embedded su( 2 ) solutionfij which are stable in 
the sphaleronic sector and fulfil (l3lT) . Happily the argument in Section (4.5) carries over identically: as before we 
use an embedded su( 2 ) solution, which is stable in both sectors given large |A|, and for which J(p) ^ 0; and we use 
analogous existence results from [T] , which we use to establish analyticity of boundary conditions and particularly 
that solutions nearby the embedded solution will also have $(p) 7 ^ 0 over the range of interest. 

At last, we deduce that if solutions stable in the sphaleronic sector exist in some neighbourhood 93 s of stable 
embedded solutions, and solutions stable in the gravitational sector exist in some neighbourhood 93 g of stable 
embedded solutions, then we can find solutions in 93s FI 9 %g which are stable in both sectors. 


4. Conclusions 

In this paper we have proven the stability of some non-trivial black hole solutions to purely magnetic topological 
su(iV) Einstein-Yang-Mills equations (with A < 0) under linear perturbations which share the same topology as the 
angular portion of the metric ©• Our focus has been to describe precisely how the change in topology affects the 
analysis carried out in [2]. 

First we applied a linear perturbation to the appropriate equilibrium field equations [lj. We found that the 
purely magnetic system decouples into two sectors. The sphaleronic sector, after some manipulation, yielded a 
hyperbolic system, requiring a positive operator for stability. We deduced that the sphaleronic sector would have 
no dynamics for the RNTadS solution, and was unconditionally stable for su(2) embedded solutions, since the 
condition for stability for su(2) topological black holes (l^Ull is guaranteed fulfilled. We then showed that genuine 
su (N) solutions could be found in some neighbourhood of these embedded solutions which have no instabilities in 
the sphaleronic sector, since the argument from [2j carries over. 

We then examined the gravitational sector. We showed that the metric perturbations could be eliminated and 
obtained a self-adjoint Schrodinger-like matrix equation in 6ojj, for which we once again required a positive definite 
matrix for stability. We based our approach on [2], using a theorem of Amann and Quittner ’s m concerning the 
number of bound states of a Schrodinger-like operator. We can again prove the existence of non-trivial solutions 
with no instabilities in the gravitational sector, in some neighbourhood of existing embedded su(2) solutions, in 
the regime of |A| —7 00 . Finally, we reasoned that there must exist some solutions in the intersection of the two 
neighbourhoods of the embedded solutions which possess no unstable modes in either sector. 

The main result of this paper is the existence of solutions to purely magnetic, topological, asymptotically 
adS su(N) field equations, in some neighbourhood of known existing solutions, which are stable under linear 
perturbations in the limit of large |A|. We also note a slightly surprising result that came from our analysis of 
the RNTadS solution. As we mentioned, we expect (by analogy with the asymptotically flat case) that for A < 0, 
at least some nodeless solutions will be stable, and therefore we suspected (by analogy with the su(2) spherical 
case [24]) that nodelessness may be a sufficient (but not necessary) condition for stability; and so in our existence 
proofs we have searched primarily for nodeless solutions as being more likely to be stable. Here we find examples of 

We note that we cannot apply the Nodal Theorem to the RNTadS solution, even though it is stable, because the gauge field in that 
case is identically zero and thus the condition on $(p) is immediately violated. 
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solutions for which nodelessness is not necessary for stability, seemingly due to the topology of the manifold. This 
discussion is outside the scope of this work but perhaps it would merit future study. 

Another possible way this work could be extended would be to use numerical methods to empirically verify the 
assertions we have made and to find some examples of stable topological solutions. Such an analysis was carried 
out for the spherically symmetric case mm, where it was verified that as |A| increased, the parameter space of 
solutions satisfying the stability condition in that case grew exponentially larger. Therefore two obvious questions 
to ask are what the space of solutions looks like in this case, and precisely how the abundance of the solutions is 
affected by the value of |A| and the N — 1 initial gauge field parameters - perhaps some bounds could be established 
connecting these quantities. 

A further natural future direction of this work would be to revisit the “No-hair” conjecture as formulated by 
Bizon [3], in light of these results, which is stated as: 

“Within a given matter model, a stable black hole is uniquely determined by global charges.” 

In light of the ‘No-hair’ conjecture, then, it can be seen that the stability of solutions is arguably more important 
an issue in the context of gravitational physics than their mere existence; but that is not the whole story, since 
there is the question of whether ‘global charges’ can be defined in this case. In [35], it is argued that spherically 
symmetric purely magnetic su(N) black holes can be characterised by a finite set of charges at infinity, and it 
would be of interest to know whether this was true for topological black holes as well, or whether the change in the 
entropy-temperature curves might cause a significant difference to the analysis. 
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